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Abstract. Infrared Space Observatory (ISO) observa- 
tions have shown that O-rich Asymptotic Giant Branch 
(AGB) stars exhibit crystalline silicate features in their 
spectra only if their mass-loss rate is higher than a cer- 
tain threshold value. Usually, this is interpreted as evi- 
dence that crystalline silicates are not present in the dust 
shells of low mass-loss rate objects. In this study, radiative 
transfer calculations have been performed to search for an 
alternative explanation to the lack of crystalline silicate 



features in the spectrum of low mass-loss rate AGB stars. 1987; van der Veen 



are referred to as Miras, whereas the stars with high mass- 
loss rates (~ 10~ 5 - ~ 10~ 3 M Q yr _1 ) are the so-called 
OH/IR stars. In the cooling outflow dust formation oc- 
curs, which adds an infrared excess to the Spectral Energy 
Distribution (SED). In case of Miras, the dust shell is op- 
tically thin, and the underlying stellar spectrum is still 
visible, but the optically thick dust shell of OH/IR stars 
completely obscures the stellar spectrum, making these 
objects only detectable at infrared wavelengths (Bcdijn 
ITabing 1988Q. 



It is shown that due to a temperature difference between 
amorphous and crystalline silicates it is possible to include 
up to 40% of crystalline silicate material in the circumstel- 
lar dust shell, without the spectra showing the characteris- 
tic spectral features. Since this implies that low mass-loss 
rate AGB stars might also form crystalline silicates and 
deposit them into the Interstellar Medium (ISM), the de- 
scribed observational selection effect may put the process 
of dust formation around AGB stars and the composition 
of the predominantly amorphous dust in the Interstellar 
Medium in a different light. Our model calculations result 
in a diagnostic tool to determine the crystallinity of an 
AGB star with a known mass-loss rate. 

Key words: stars: AGB and post-AGB - circumstel- 
lar matter - dust, extinction - infrared: stars - radiative 
transfer 



1. Introduction 



Important reasons for studying the nature of dust in 
AGB stars can be found in the fact that these stars are 
the main contributors of dust into the interstellar medium. 
The properties of the dust in the ISM can thus be directly 
compared to the dust properties in shells around evolved 
stars. This will provide information on the processing of 
the dust in the ISM and will improve our insight into the 
life cycle of dust. Moreover, studying the condensation of 
dust in the outflows of AGB stars will provide useful diag- 
nostics for the determination of the physical conditions in 
these outflows, thus increasing our understanding of this 
important stage of stellar evolution. 

Most AGB stars exhibit an oxygen-rich chemistry, 
i.e. a number ratio C/O < 1, in their dust shells, which 
gives rise to the condensation of oxygen-rich dust species. 
Based on early infrared observations, the presence of sil- 
icate species has been proposed by several authors (e.g. 
pilman 1969} |Joncs fc Merrill 1976| ). Pre-ISO studies al- 
most exclusively pointed toward the presence of amor- 
phous Fe-bearing silicates; such dust shells are for exam- 
plc studied by [Butchart fc Whittet (1983| ) 



Pegourie & 



Papoular (1985| ) and Bcdijn (1987). Infrared spectroscopy 



Oxygcf-rich Asymptotic Giant Branch (AGB) stars are 

M-tvn e stars of intermediate mass (1 M c < M < RM C A in performed w ith the Infrared Space Observatory flKesslei 



the lat k stages of stellar evolution which exhibit He-' and ct g jjjg§ revealed the presence of crystalline silicates 



H-shell burning (for a review on AGB stars, see Habing 
|(1996| ), and references herein). These stars undergo an in 
creasing mass loss while they evolve along the AGB; the 
stars with low mass- loss rates (~ 10 -8 - ~ 10~ 6 M Q yr _1 ) 



next to the amorphous silicate dust in the circumstcllar 
dust shells of some evolved stars (Waters et al. 1996). In 
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addition, the observed wavelengths of the emission bands 
of the crystalline silicates in evolved stars suggest that 
these grains contain little or no Fe, but are very Mg-rich 
(Molstcr ct al. 1999). The data are consistent with no Fe 
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in the grains. There seems to be a correlation between 
the peak of the SED, i.e. the optical depth toward the 
central star, and the presence of spectral features of crys- 
talline olivine (forsterite, Mg2SiO/t) and crystalline pyrox- 
ene (enstatitc, MgSiC^) (|Waters et al. 1996|; |Cami et al 



spherical grains, for which we use Mie calculations to de- 
termine the absorption and scattering coefficients. The 
adopted chemical composition and laboratory measure- 
ments of optical constants - necessary to calculate the ex- 
will be discussed in Sect 



tinction properties 



2.1 



1998 |rpylvcstcr ct al. 1999| ), in the sense that this material size distribution we used a powerlaw N(a) oc a 3 5 , con 



For the 

3.5 



is not seen in stars that are not obscured. This is usually sistcnt with the grain-size distribution of the ISM (Mathis 



interpreted as evidence that crystalline silicates are not et al. 1977). Following studies on the sizes of interstellar 



presen t in the dust shells around low mass-loss rate AGB grains and of S t ardust recovere d from meteorites (Andcr, 



stars, which have an optically thin dust shell. Apparently, 
a threshold value for the mass-loss rate M is required to 
enable the condensation of crystalline silicate in the stellar 
outflows. This critical density would thus put constraints 
on the dust condensation sequence. Various explanations 



hav e been offered for this threshold density ( Fielens et al 
1998 |; |(j!ail fc Sedlmayr 1999| ; |Sogawa fc Kozasa 1999| ). Yet 



its existence is not well established. 

In this study, we will evaluate the existence of the 
threshold value for the mass-loss rate below which crys- 
talline silicates are not detected. Does the lack of crys- 
talline silicate features in the SED of Miras really imply 
that the crystalline silicates are not present, or is some 
other effect at work suppressing the characteristic narrow 
features? We will study this issue by calculating a grid of 
detailed model spectra of circumstellar dust shells for a 
wide range of mass loss and degree of crystallinity. 

In Sect. El the model for the circumstellar dust shell 



fc Zinner 1993| ; Kim et al. 1994| ) we have adopted a grain- 
size range of 0.01 - 1 /im. The dust shell consists of sep- 
arate grain populations, i.e. each dust grain contains one 
type of silicate, implying that the separate dust species 
can have different temperature profiles. Finally, all dust 
species have the same spatial distribution. 

The models in this study cover a wide range in mass- 
loss rates and degree of crystallinity. The mass-loss rate is 
varied from 5 • 10~ 8 to 10~ 4 Mq yr _1 . The crystallinity x, 
defined as 

total mass crystalline silicates 
total mass silicates 

is varied from 0% to 50%. The dust shell is illuminated at 
the inner radius by a typical M9 giant, of which the spec- 
trum is taken from Fluks et al. (1994). The size of the inner 
radius i?; n of the dust shell is determined by the conden- 
sation temperature of the dust, T C ond ~ 1000 K (Gail & 



will b d described, including the applied optical constant 
for the silicates. In Sect. |^ the results will be discussed for 
varying mass-loss rate and degree of crystallinity. Sect, [l] 
contains the physical explanation of the obtained results 
and highlights possible implications for the evolution of 
dust. A summary of this study is given in Sect. |[ 

2. Modelling the circumstellar dust shell 

The dust radiative transfer programme modust (de Koter 
et al., in prep.) has been used to calculate detailed spec- 
tra of circumstellar dust shells. The code solves the trans- 
fer equations subject to the constraint of radiative equi- 
librium, fixing the temperature distribution of the dust 
grains. The solution technique focusses on the combined 
moment form of the transfer equation, such that the scat- 
tering term in the source function can be solved for ex- 
plicitely, updating the Eddington factors throughout the 
lambda iteration by means of a ray-by-ray formal solution. 

In the mode used in this paper, the dust is assumed 
to be distributed in a spherical shell with a radial density 
profile p cx r~ 2 , i.e. that of a time-independent stellar 
wind outflow at constant velocity 

M = 47T7- 2 p{r)v cxp 

where M is the mass-loss rate and v exp the outflow veloc- 
ity. The program allows one to specify an arbitrary num- 
ber of dust specie s, each with its own s hape and size dis- 
tribution (see also Bouwman et al. 200C ) . Here we adopted 



Sedlmayr 1999]) . The outer radius is set to i? Q ut = 200i? 



which is sufficiently far out to account for all the near in- 
frared (NIR) and mid infrared (MIR) flux. The outflow 
velocity is assumed to be at a constant value of v, 



cxp 



20 

km s _1 . The dust/gas mass ratio in the outflow is taken 
to be / = 0.01. The last two numbers are typical for Miras 
and OH/IR stars. The presented results are expressed in 
terms of mass-loss rates for the above value of the outflow 
velocity. However the critical parameter determining the 
spectrum is obviously the dust density pdust(^)- Therefore 
the results are invariant for the quantity 



Q = 



M 



& for /«! 



which follows from the mass continuity equation. 



2.1. Dust optical constants 

Previous studies of the SED of dust shells surrounding 
AGB stars Qedijn 198~7|; [Justtanont fc Ticlcns 1992| ; g 



|Si(lancr fc Lc Bcrtre 1993 , 199E) used so-called astronom- 
ical or dirty silicates, of which the optical constants were 



observationally determined (Jones & Merrill 1976; Drainc 
|fc Lee 1984 ). The chemical composition of astronomical 
silicate is not known, although likely the extra absorption 
is due to the presence of iron in the form of Fe 2+ in the 



mineral (Schutte fc Tielens 1989) 



In the last few years, accurate laboratory measure- 
ments of the optical constants of amorphous silicates 
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(Dorschner et al. 1995) and crystalline silicates flKoike 
et al. 1 |993| ; |Koike fc Shibai 199S| ; |Jager et al. 1998| ) have 
become available. Together with the high spectral reso- 
l ution of the ISO Short Wavelength Spectrometer (SWS) 
dde"& aauw et al. 1996|) and L ong Wavelength Spectrome- 
ter (LWS) flClegg et al. 1996| ) this opens the possibility to 
determine the dust composition in great detail, at least in 
principle. However, model calculations have shown that 
the SED can not be fitted using only the optical con- 
stants of amorphous and crystalline silicates. An addi- 
tional source of opacity in the NIR appears to be miss- 
ing (Kemper et al., in pr ep). This is a well-known prob- 
lem which dates back to Jones fe Merrill (1976 ); Bcdijij 
(19870 and |Schuttc fc Ticlcns (1989| ). Although it is a 



fundamental issue, for the present study it suffices to use 
the optical constants of the amorphous sili cates calc ulated 
from the observed spectra of AGB stars ( Suh 199S ). The 
thus obtained optical constants of the AGB silicate can- 
not provide information on the exact composition of the 
amorphous component of the dust shell, however for our 
studies of the appearance and strength of the crystalline 
silicate features, the use of the "Suh-silicate" will not af- 
fect our results. As the peak position of the crystalline 
features in laboratory samples is strongly dependent on 
the Fe-content, ISO observations of AGB and post-AGB 
stars demonstrate that the dust shells only contain Mg- 



rich crystalline silicates (Molster 
et al. 



et al. 1999; Sylvester 



l|999|). For these crystalline silicates, we used optical 
constants of two different types: forsterite (Mg2Si04) and 



enstatite (MgSiOa) (Jager et al. 1998). It is important to 



realize that because of the low Fe-content, their absorp- 
tivity in the Near Infrared (NIR) is very low. Equal mass 
fractions of both crystalline species have been used, unless 
stated otherwise. 



3. Model results 

Here we present the model results for increasing mass- 
loss rates (Sect. 3.1) and increasing degree of crystallinity 
(Sect. [3~2] ). In Sect. |3.3| the observational constraints of 
ISO will be taken into account and applied to the model 
spectra. 



3.1. Varying the mass-loss rate 

In Figs. [j], H and |^ results of the model calculations are pre- 
sented. Each figure shows the fluxes F„(A[/im]) for mass- 
loss rates varying from 10~ 7 to 10~ 4 Mq yr _1 for a con- 
stant degree of crystallinity. In Fig. [l] results for a dust 
shell consisting of completely amorphous dust are shown. 
Clearly visible are the bands due to Si-0 stretching and 
O-Si-O-bending around 10 and 20 /im respectively. The 
strength of these bands varies with increasing mass-loss 
rate; for low M the dust shell is optically thin and the 
amorphous silicate bands appear in emission. The strength 
of the emission bands increases with increasing mass-loss 
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Fig. 4. Contrast plot of the 25 - 50 /im region for dif- 
ferent mass-loss rates at a constant degree of crystallinity 
x = 0.10. The M (Mq yr -1 ) values are indicated for each 
curve. The contrast is obtained by dividing the model 
spectrum containing 10% crystalline dust by the model 
spectrum of a dust shell with an equal amount of amor- 
phous dust, but lacking the crystalline dust component. 

rate. When the radial optical depth of the peaks of the 
spectral features approaches unity, the features become 
self-absorbed. This occurs around M = 10~ 5 Mq yr _1 
in case of the 10 fim feature. For the 20 /im feature this 
transition is less well defined, and occurs at a somewhat 
higher mass-loss rate. 

Figs. |^ and |3| show the emerging spectra for a crys- 
tallinity of 10% and 20% respectively. The overall shape 
of these spectra resemble that of the amorphous silicate 
dust shell spectra in Fig. [I]. However, while the degree of 
crystallinity remains constant, the narrow crystalline sil- 
icate features at 23.6, 27.6, 32.5, 36.1, 41.2 and 43.2 /im 
become quite apparent at relatively high mass-loss rates. 

In the models with low mass-loss rate, the narrow crys- 
talline silicate features are not discernible. Let us first in- 
vestigate the relative strength of the crystalline silicate 
features with respect to the amorphous dust continuum 
in some more detail. For this purpose the contrast is cal- 
culated by dividing the spectrum of a dust model of a 
given crystallinity by the spectrum of a model with equal 
properties and the same amorphous silicate dust mass, 
but which is lacking the mass component taken up by the 
crystalline silicates. Fig. ^ shows the contrast for three dif- 
ferent mass-loss rates, including 10% crystalline material. 
The narrow bands due to crystalline silicates at 27.6, 32.5, 
36.1, 41.2 and 43.2 /im show that the relative strength in- 
creases with mass-loss rate. Note that the continuum be- 
tween the crystalline silicate features does not necessarily 
return to unity, as the "extra" IR emission from crystalline 
silicates raises the temperature of the amorphous silicates 
slightly. For M = 10~ 7 Mq yr _1 the contrast between the 
crystalline silicate peaks and the surrounding continuum 
is less than 1% for all features; for M = 10~ 6 Mq yr" 1 
some peaks have a relative strength of ~ 4% above contin- 



4 



F. Kemper et al.: Crystallinity versus mass- loss rate in Asymptotic Giant Branch stars 




100 



Fig. 1. Fluxes calculated for a 
completely amorphous dust shell, 
i.e. x — 0, for a wide range of mass- 
loss rates, calculated for v cxp = 20 
km s _1 and / = 0.01. In each panel, 
the mass-loss rate M (Mq yr" 1 ) 
is indicated. The spectra are nor- 
malised on the maximum value in 
the 2 - 100 /im wavelength region. 
Around 10 and 20 fxm broad fea- 
tures due to amorphous silicates are 
present in the spectra, occurring ei- 
ther in emission or absorption, de- 
pending on M (see text). Note that 
in case of optically thin dust shells 
(low M) the stellar spectrum with 
its characteristic molecular absorp- 
tion bands is still visible at A < 8 
yitm, whereas for the optically thick 
case (high M), the stellar flux at 
NIR wavelengths is completely ab- 
sorbed by the dust. 




uum; whereas the strongest peaks of the M = 10 -5 M & 
yr -1 model show a contrast of ~ 10% above continuum. 
The detection limit of ISO for broad spectral features in 
this wavelength region is — depending on the quality of 
the spectrum — around 5% of the continuum level, which 



Fig. 2. Normalised fluxes calculated 
for a dust shell containing silicates 
with a degree of crystallinity x = 
0.10, for the same mass-loss rates as 
in Fig. 0. The narrow features due to 
crystalline silicates start to appear 
in the 20 - 50 jxm range are only 
appearing for relatively high mass- 
loss rates. At M » 10" 5 M Q yr" 1 
the crystalline silicates are also dis- 
cernible as substructure in the 10 
/j,m feature (see also Fig. |ll]) 

implies that of the three curves plotted in Fig. only the 
ISO spectrum of a star with an outflow of 10~°Mg yr -1 
would result in a unambiguous detection of crystalline sili- 
cate features. A detection in the spectrum of an AGB star 
with M = 10~ 6 Mq yr" 1 would be dubious. A typical 
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100 



Fig. 3. Normalised fluxes calculated 
for a dust shell containing silicates 
with a degree of crystallinity x = 
0.20, for the same mass-loss rates as 
in Fig. [l]. The narrow features due 
to crystalline silicates in the 20 - 
50 //m range are only appearing for 
relatively high mass-loss rates. At 
M w 5 • 10" 6 Mq yr- 1 the crys- 
talline silicates are also discernible 
as substructure in the 10 /im feature 
(see also Fig. Ill]) 



Mira, with M = 10~ 7 M© yr -1 , with a crystallinity of 
10% would not show the characteristic features. This re- 
sult is consistent with the threshold mass-loss rate above 
which crystalline silicate s are detected with ISO observa- 
tions ( Waters ct al. 199(j ; |Cami ct al. 1998 ; Sylvester et al 
19990 



3. 2. Varying the degree of crystallinity 

In order to determine the amount of crystalline silicates 
that should be present in a typical Mira dust shell such 
that the narrow features would be detected in an ISO 
spectrum, model spectra with different degrees of crys- 
tallinity have been calculated. For a typical Mira mass- 
loss of M — 10~ 7 Mq yr -1 the crystallinity has been 
varied from 5% to 50%. The spectra of the x — 0.10 and 
x = 0.20 models, discussed in Sect. 3.1, are included in the 



overview. The contrast has been calculated in the same 
way as described above. The result is shown in Fig. S, 
for models with M — 10~ 7 Mq yr -1 . The dashed line 
indicates the 5% detection limit of ISO for solid state fea- 
tures, typical for A ~ 25 - 50 /im, and thus gives an im- 
pression of which features will be discernible. As stated 
before, the relative flux of the crystalline silicate model 
in between the narrow crystalline silicate features is not 
necessarily equal to that of the amorphous model fluxes. 
The peak strength should therefore be compared to the 
local continuum, rather than to unity. For this purpose an 
additional base-line has to be subtracted. Following this 
method one easily sees that very high mass abundances of 




30 40 
Wavelength [jttm] 



50 60 



Fig. 5. Contrast plot for A < 60 /im for different degrees 
of crystallinity x at a constant mass-loss rate of M = 10~ 7 
Mq yr -1 . The crystalline fraction x is indicated at each 
curve. The contrast is obtained by dividing the model 
spectrum of a dust shell with a certain degree of crys- 
tallinity by the model spectrum of a dust shell with ex- 
actly the same amount of amorphous dust, but which lacks 
the crystalline dust component. The ISO detection limit 
for emission features of 5% with respect to a continuum at 
^cryst+amorph/^amorph = 1 is indicated with a dashed 
line. 



crystalline silicates are required for a marginal detection 
of crystalline silicate features in a typical Mira spectrum. 
For a crystallinity of x = 0.40 only the 43.2 /im and 32.5 
/xm features can be detected by ISO. 
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3.3. A trend in crystallinity? 

From Sects. [O] and 3.2 it becomes clear that the com- 
parison of the ISO detection limit with the contrast of 
the crystalline silicate features is essential for studying 
the correlation between mass-loss rate and crystallinity 
of AGB stars. From Fig. || the 43.2 /jm feature due to 
crystalline pyroxene (enstatite, MgSiOa) has been selected 
for closer examination, which is a conservative choice. It 
shows the largest contrast with the local continuum, and 
would therefore be the first feature to be detected in the 
spectrum when mass loss or crystallinity increases, under 
the assumption that the mass fractions of crystalline py- 
roxenci and foroterito are equal. 



Table 1. Relative fluxes and mass- loss rates of the sample 
stars from Sylvester et al. (1999| ) 



Star 


rel. flux 
32.5 (im 


rel. flux 
43.2 /im 


M 
M @ yr" 


i 


ref. 


o Cet 


< 0.05 


< 0.05 


1.0 


■ io- 


7 


1 


CRL 2199 


< 0.05 


< 0.05 


1.5 


■ 10' 


5 


2.3 


WX Psc 


0.115 


0.090 


1.9 


■ 10" 


5 


2,4,5,6 


OH 104.9+2.4 


0.064 


< 0.05 


6.0 


■ 10" 


5 


6 


OH 127.8+0.0 


< 0.05 


0.15 


2.0 


■ io- 


4 


3 


OH 26.5+0.6 


0.062 


0.095 


1.6 


■ io- 


4 


3,6,7 


AFGL 5379 


0.077 


0.093 










OH 32.8-0.3 


< 0.05 


0.16 


1.0 


■ 10" 


4 


3,6,8 



1 [Planesas et al ( 19901) 2 |Knapp 



■ fe Morris 



1985!) 3 Uusttanont 



Tielens (19921) 4|Sopka et al. (198% 5 te Bertre fe Winters 



For a 1'allgc of MCTM MB khd M gjj| 6 | s ^ utt ' c fc ' gj^ (i 98S | ) 7 frusttanont et aL7T99tf 



rates the contrast spectra have been calculated. To de- 
termine the relative peak strength with respect to the lo- 
cal continuum, a first order polynomial base-line has been 
subtracted from the feature, such that the local contin- 
uum is equal to unity. Then the peak value of the 43.2 /im 
feature was measured, and plotted in Fig. ||. Note that the 
enstatite mass fractions differ a factor 2 with the elsewhere 
used crystallinity x, because it was assumed in Sect, p.l 
that enstatite and forsterite are equally abundant. Fig. p 
is a diagnostic tool to determine the crystalline pyroxene 
content of the silicates present in the dust shell around an 
AGB star. To do so, one has to determine the relative peak 
strength of the 43.2 /xm feature with respect to the local 
continuum, a quantity plotted on the vertical axis, and 
combine this with the mass-loss rate of the object. The 
data points of the sample of AGB stars previously anal- 
ysed by Sylvester et al. (1999] ) are presented to illustrate 
the method. The relative strength is measured from ISO 
spectra with an accuracy of 2 per cent of the continuum 
(Table 0). We adopted an accuracy of a factor 2 for the 
mass-loss rates. Using these data, one can estimate from 
Fig. U that the enstatite makes up around 3% ± 2% of the 
total silicate dust mass around WX Psc and OH 26.5+0.6, 
and around 6% ± 2% in OH 127.8+0.0 and OH 32.8-0.3. 
In CRL 2199 and OH 104.9+2.4 the enstatite mass frac- 
tion is well below 2-3%, whereas the upper limit obtained 
for o Cet (Mira itself) does not put interesting constraints 
on the enstatite content. The dust shell of o Cet could 
still easily contain a mass fraction of 20% of crystalline 
pyroxene. A constant enstatite mass fraction of ~ 5% is 
consistent with the observed strength of the 43.2 /im fea- 
ture for a wide range of mass-loss rates. 

A similar analysis can be performed on the forsterite 
(Mg2Si04, crystalline olivine) mass fraction, using the 
32.5 /xm feature, which has the best contrast. For the 
same range of mass-loss rates the emerging dust spectrum 
has been calculated. The mass fraction of the crystalline 
olivine varied from 0% to 20%. All other input parame- 
ters had equal values as the parameters in the analysis of 
the crystalline pyroxene mass fraction. We calculated the 
relative peak value of the 32.5 /im feature, by subtract- 



8 GIroenewegen (1994) 



0.60 



0.50 



a. 



0.40 r 



0.30 r 



t 0.20 t 



0.10 



0.00 




10 ' 10 ' J 10 

mass-loss rate (M a yr -1 ) 

Fig. 6. Relative strength of the 43.2 /im feature with re- 
spect to the local continuum as determined by amorphous 
material as a function of M. The curves are calculated for 
enstatite mass fractions of 0.05; 0.10; 0.15 and 0.20; which 
are indicated in the plot. The dashed line represents the 
ISO detection limit at 43.2 /im, which is ~ 5%. For several 



AGB stars, M and F, 



ZAcryst 



are indicated (see Ta- 



ble |l|). For some stars, only an upper limit to the relative 
flux is given. 



ing a first order polynomial from the forsterite contrast 
plots. Fig. shows the relative peak fluxes as a function 
of mass-loss rate, calculated for different mass fractions of 
forsterite. From this plot it is obvious that in low mass- 
loss rate AGB stars, the strength of the 32.5 /im feature 
is less than 5% of the continuum, which makes it imper- 
ceptible for ISO, as is illustrated with the data for o Cet. 
For higher mass- loss rate AGB stars, the strength of the 
feature provides a reliable measure to determine the crys- 
talline olivine mass fraction. The forsterite mass fraction 
of our sample of high mass-loss stars however, is at most 
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Fig. 7. Relative strength of the 32.5 /im forsterite feature 
with respect to the local continuum as a function of M. 
The curves are calculated for forsterite mass fractions of 
0.05; 0.10; 0.15 and 0.20; which are indicated in the plot. 
The dashed line represents the ISO detection limit at 32.5 
/tm, which is ~ 5%. For a sample of AGB stars, M and 
Fv,ayvt/Fu,cant are indicated (see Table 0). For some stars, 
only an upper limit to the relative strength is available. 



5% ± 2% for WX Psc. There is no evidence for an increas- 
ing forsterite mass fraction for higher mass-loss rates. 

Figs. H and [)] illustrate that an observational selection 
effect is present in the correlation between mass-loss rate 



and crystallinity reported in previous studies by Waters 
et al. Q1996D ; |Cami et al. (1998| ); |Sylvcster et al. (1999] ) 
The OH/IR stars with M > 10~ 5 M© yr" 1 only have to 
contain a small mass fraction of crystalline pyroxene in 
their dust shells in order to exhibit the 43.2 /im or 32.5 
/tm feature. A few percent is already perceptible with ISO. 
However, the 43.2 and 32.5 /tm features in the spectra of 
Miras, with M < 10~ 6 M Q yr" 1 , will only be detected 
for high crystallinity. For example, in the dust shell of o 
Cet, with M w 10~ 7 Af Q yr" 1 , a dust mass fraction of 
20% of enstatite could be present without showing the 43 
fim feature in its spectrum. The same goes for the 32.5 
/im forsterite feature. If we assume that equal amounts 
of crystalline pyroxene and olivine are present, the ISO 
detection limit only provides an upper-limit to the crys- 
tallinity of x = 0.40. If a constant crystallinity of x = 0.10 
(i.e. a mass fraction for both enstatite and forsterite of 
5% of the total silicate mass) is assumed, while the star 
evolves along the AGB, the 32.5 fim and 43.2 /mi fea- 
tures will appear in the spectrum when the mass loss is 
sufficiently high. The threshold value, defined by the ISO 
detection limit, is found at M ~ 10~ 6 M yr" 1 . This is 
consistent with observational studies that conclude that 
crystalline silicate features only appear to be present in 
high mass-loss rate objects. 



Fig. 8. Relative peak strength of the 43.2 /im feature as 
a function of mass-loss rate and fraction of mass taken by 
enstatite (MgSiOs). The peak strength is determined with 
respect to the local continuum. A relative peak strength 
> 5% is observable with ISO. The shaded area represents 
the values for M and enstatite mass fraction for which the 
strength of the 43 /tm feature is below the ISO detection 
limit. The mass fraction of enstatite is sampled at 0.00, 
0.05, 0.10 and 0.20; and the mass-loss rate is sampled at 
M/10 6 = 0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50 and 100 M Q 
yr _1 . 
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Fig. 9. Relative peak strength of the 32.5 /tm feature as 
a function of mass-loss rate and fraction of mass taken 
by forsterite (Mg2Si04, crystalline olivine). The peak 
strength is determined with respect to the local contin- 
uum. A relative peak strength > 5% is observable with 
ISO. The shaded area represents the values for M and 
forsterite mass fraction, for which the strength of the 32.5 
/tm feature is below the ISO detection limit. The mass 
fraction of forsterite is sampled at 0.00, 0.05, 0.10, 0.15 
and 0.20; and the mass-loss rate is sampled at M /10 6 = 
0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50 and 100 M yr" 1 . 
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4. Physical explanation and discussion 

The appearance of crystalline silicate features with in- 
creasing mass-loss rate for a constant degree of crys- 
tallinity may be caused by the difference in absorptivity 
in the NIR of crystalline and amorphous silicates. This 
difference is due to a difference in Fe-content of the amor- 
phous and crystalline components. The crystalline sili- 
cates investigated do not contain any iron, whereas in the 
amorphous silicates the iron and magnesium contents are 



roughly equal. Previous laboratory studies (Koike et al 



1993; lager et al. 1998) have shown that the peak po- 



sitions of crystalline silicate features strongly depend on 
the Fe-content of the silicate. At present, there have not 
been any detections of crystalline silicates containing iron 
around evolved stars, which justifies the use of laboratory 
data of pure Mg 2 Si04 and MgSiC>3. 

The absorptivity in the NIR and at shorter wave- 
lengths increases with increasing Fe-content. The amor- 
phous silicates thus absorb more radiation from the central 
star, which emits most of its flux at NIR wavelengths, than 
the crystalline silicates do, resulting in a large tempera- 
ture difference between amorphous and crystalline mate- 
rial, the latter staying much cooler. In the mid and far 
infrared the absorptivity of both dust types are more or 
less equal. In an optically thick dust shell, where the stellar 
flux is absorbed and re-radiated several times by the dust, 
the grains in the outer regions are illuminated by a radia- 
tion field that typically peaks at mid infrared wavelengths. 
As in this wavelength regime only a small difference in ab- 
sorptivity occurs, a large temperature difference between 
the amorphous and crystalline dust components will not 
develop. Summarising, at the inner radius of an optically 
thick dust shell the crystalline silicates are much colder 
than the amorphous silicate grains, but the temperature 
difference decreases at larger distances. 

The implications for the strength of the crystalline sil- 
icate features due to this temperature difference are quite 
dramatic. In low mass- loss rate AGB stars (Miras), the 
dust shell is optically thin, implying the entire dust shell 
is visible and contributes to the SED. The relatively hot 
amorphous silicates at the inner radius then dominate the 
SED of a Mira, and the contrast of the crystalline silicates 
with respect to the amorphous silicate spectrum becomes 
very poor. In case of a high mass-loss rate AGB star, only 
the outer regions of the dust shell are visible, because the 
shell is optically thick at IR wavelengths. The tempera- 
ture differences in the outer layers are small, as explained 
above, and therefore the crystalline and amorphous dust 
components both contribute in comparable amounts to 
the spectrum. This improved contrast of the crystalline 
silicates with respect to the amorphous dust emission en- 
ables the detection of the narrow crystalline features in 
the spectrum of OH/IR stars. 

The ISO detection limit allows the presence of a sig- 
nificant fraction of crystalline silicates in the circumstel- 
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Fig. 10. ISO SWS spectrum of AFGL 2999. The spectrum 
of this star is typical for that of an object somewhere in 
between typical Miras and OH/IR stars. High resolution 
and high signal-to-noise observations may perhaps reveal 
crystalline silicate features as sub-structure in the broad 
amorphous 10 /im feature. No crystalline silicates are seen 
at the 5 percent level in the 20 - 45 /im range. The artefact 
at 28 /im is due to calibration problems with band 3e of 
ISO SWS. 



lar shells of low mass-loss rate AGB stars. In fact, the ISO 
data are consistent with model calculations assuming that 
the crystallinity of the dust is constant with mass- loss. 



Until now, most theoretical studies (Ticlcns et al. 1998 
|Gail fc Scdlmayr 199S ; Sogawa fc Kozasa 1999] ) assumed 



that crystalline silicates are only formed in high density 
outflows. The results presented in this study suggest that 
the possibility of the condensation of crystalline silicates 
in low mass-loss rate AGB stars should be considered in 
future studies. 

4-.1. How to search for crystalline silicates in Miras 

To date, there are no significant detections of crystalline 
silicates in low mass-loss rate AGB stars reported. As 
shown above, this provides a relatively high upper limit to 
the mass fraction taken by crystalline silicates in the shells 
around these stars. In order to lower this upper limit and 
to investigate the presence of crystalline silicates in Miras, 
additional observations with a better detection limit than 
that of ISO, are required. Unfortunately, the most promi- 
nent and characteristic features of crystalline silicates ap- 
pear in the 25 - 50 /^m region, and with the end of the ISO 
mission the opportunities to observe in this wavelength 
region are limited. An important future mission covering 
part of this wavelength region is SIRTF, which can per- 
form intermediate resolution spectroscopy (R ~ 600) up 
to A « 37 /im. The crystalline silicate feature with the 
largest contrast value within this wavelength range is the 
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Fig. 11. Predicted flux levels in the 7-13 /im region for 
an AGB star with M = 5 ■ 10~ 6 M© yr -1 . Model spec- 
tra with three different crystallinities (x — 0.0; x — 0.10; 
x = 0.20) are indicated. The broad emission feature is due 
to amorphous silicates. The substructure is due to crys- 
talline silicate features that appear in absorption. The ISO 
SWS spectrum of AFGL 2999 is plotted for reference; in 
this spectrum it is not possible to detect crystalline sili- 
cate features, since the the substructure in the spectrum 
is mostly due to observational noise. 



32.5 /im forsterite feature. In Fig. ||, the contrast of this 
feature with respect to the continuum is presented as a 
function of mass-loss rate and Mg2SiC>4 mass fraction. 

In addition, it is also possible to study the 10 and 
20 /im silicate features, which both can be observed from 
the ground. Since it is already known that the high mass- 
loss rate OH/IR stars contain crystalline silicates in their 
circumstellar dust shells, it is particularly interesting to 
study the 10 /im feature of low and intermediate mass- 
loss rate objects, with M < 10~ 5 Mq yr -1 . From figs. || 
and |^ it one can conclude that the 10 /im feature of in- 
termediate mass- loss rate AGB stars (M « 10~ 6 to 10~ 5 
Mq yr -1 ) would be the most likely to reveal the presence 
of crystalline silicates. The 10 /im region has the intrin- 
sically strongest resonance peaks of crystalline silicates. 
Superposed on the 10 /im feature of amorphous silicates 
which is just optically thick (r > 1) these strong reso- 
nances appear as absorption features in the spectrum. The 
intermediate mass-loss rate AGB stars can be easily recog- 
nised from the ratio between the strength of the 10 and 20 
/im features. For example AFGL 2999, of which the ISO 
SWS spectrum is shown in Fig. [l0|, can be identified as a 
M = 5 • 10~ 6 M Q yr" 1 AGB star. 

Fig. p"T| shows predicted 10 /im features. The flux levels 
are calculated for models with M = 5 ■ 10~ 6 Mq yr -1 and 
x — 0.00, 0.10, 0.20. The broad amorphous silicate feature 
is partially optically thick and shows a flattened profile 
due to self-absorption. The additional opacity provided 



by the crystalline silicates appears as narrow absorption 
superposed on the broad amorphous silicate feature. The 
ISO spectrum of AFGL 2999 is plotted for reference, but 
the substructure in this spectrum is observational noise 
and it is therefore hard to disentangle the crystalline sili- 
cate mass fraction. However, it would certainly be possible 
reach a high S/N ratio and detect a crystallinity of 10% 
with the current ground-based spectrometers that cover 
the 10 /im region, such as TIMMI2 mounted on the 3.6m 
ESO telescope. 



4-2. Implications for the dust formation and processing 
in the ISM 

AGB stars are considered the main contributors of dust to 
the ISM, so one would expect a match between the dust 
composition in the ISM and in the shells of AGB stars. 
However, despite all evidence for crystalline silicates in 
AGB stars, there has as yet not been found any reliable 
evidence for crystalline silicates in the line of sight toward 
the Galactic Center (GC) ( Lutz et al. 1996 ), or other lines 
of sight through the ISM ( [Demyk et al. 1999] ). In fact, the 
line of sight toward the GC is found to resemble closely 
the silicate feature of the evolved star /i Cep, and can be 



fitted with completely amorphous olivine (Vricnd 1999, 
Bouwman et al. in prep). Cesarsky et al. (2000) report a 



questionable detection of crystalline silicates in the Orion 
Bar. 

On the other hand, it is doubtful that Miras contribute 
significantly to the dust population in the ISM. It is gen- 
erally accepted that the Mira phase lasts ~ 3 • 10 5 yr, and 
the OH/IR phase - 3000 yr flHabing 1996| ). The mass- 
loss rates for Miras and OH/IR stars (AfMira ~ 10 -7 Mq 
yr -1 , and M r/ir ~ 10~ 4 Mq yr -1 ) can be used to es- 
timate the total mass lost in each of the two phases; for 
OH/IR stars this is ~ 0.3Af Q , for Miras ~ 0.03M Q . Un- 
der the assumption that all Miras will eventually evolve 
into OH/IR stars, only a minor fraction (~ 10%) of the 
interstellar dust population is formed in a Mira dust shell. 
So, the possibility that Miras contain a significant frac- 
tion of crystalline silicates in their dust shells does not 
heavily aggravate the discrepancy between the dust com- 
position in the ISM and that in AGB dust shells, under the 
assumption that all Miras eventually evolve into OH/IR 
stars. However, if Mira stars indeed deposit dust with a 
high degree of crystallinity into the ISM, and a signifi- 
cant fraction of Miras does not evolve into OH/IR stars, 
it is more surprising that the interstellar silicate grains 
appear to have an amorphous lattice structure. Grain- 
grain collisions and ion bombardments might cause the re- 
quired amorphitization. Due to the collision, the grain will 
partially melt and consequently solidify in an amorphous 
f orm on timcscales sh orter than the annealing timescales 
( Sylvester et al. 1999 ). In a future study (Kemper et al., 
in prep.) we will thoroughly examine the differences be- 
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tween silicates in the ISM and in circumstellar shells, and 
discuss possible mechanisms of grain processing. 

5. Summary 

In this work, three important aspects of crystalline sili- 
cates around AGB stars have been studied by modelling 
the dust shell. First, we have calculated spectra emerg- 
ing from circumstellar dust shells for different mass-loss 
rates. Moreover, the effects of different degrees of crys- 
tallinity on the SED have been taken into account, as a 
function of optical depth toward the central star. Finally, 
we provide a useful and easy to use diagnostic tool to de- 
termine the degree of crystallinity for an AGB star with 
a known mass-loss rate. 

The dust shells of low mass-loss rate AGB stars can 
contain a significant fraction of crystalline silicates, while 
the characteristic sharp spectral peaks are indiscernible 
with ISO. This is caused by a large temperature differ- 
ence in the inner parts of these optically thin dust shells, 
such that the relatively warm amorphous silicate emission 
dominates the SED. The temperature difference between 
the warm amorphous dust and the cold crystalline dust is 
caused by a difference in absorptivity at NIR and visible 
wavelengths, where the amorphous dust component more 
efficiently absorbs the stellar radiation. The low absorp- 
tivity of the crystalline component is due to the very low 
Fe-content of these dust species. In case of high mass-loss 
AGB stars, the radiation from the central star is absorbed 
by the dust and re-radiated several times, so that the dust 
in the outer parts receives a relatively red radiation field, 
which peaks in the mid-IR. In that region the absorptivity 
of crystalline and amorphous silicates are similar, leading 
to similar temperatures of both species. The dust shell is 
optically thick, so only the outer layers, where crystalline 
and amorphous silicates have similar temperatures, are 
visible, leading to a high contrast of the crystalline sili- 
cate features. 

From these results one may conclude that crystallinity 
is not necessarily a function of mass-loss rate, as many 
observational and theoretical studies have suggested. The 
threshold value for the mass-loss rate above which crys- 
talline silicates are observed is consistent with the ISO 
detection limit if a constant crystallinity of 10% is as- 
sumed. The strong observational selection effect under- 
mines previously drawn conclusions regarding this thresh- 
old value. The meaning and the existence of the threshold 
mass-loss rate should be reconsidered. Future observations 
of the 10 /zm silicate feature of low mass-loss AGB stars 
at high spectral resolution can probably further constrain 
the upper limit for the crystalline silicate content in Miras. 
The low degree of crystallinity of grains in the interstellar 
medium with respect to grains in AGB shells, however, 
remains subject to further study. 
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